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ABSTRACT 

■ We report on the amazing y-ray activity from the flat spectrum radio quasar (FSRQ) PKS 1510-089 observed by the AGILE satellite in March 

1 I 2009. In the same period a radio-to-optical monitoring of the source was provided by the GASP-WEBT and REM facilities. In the radio band we 

made use also of multi-epoch 15-GHz Very Long Baseline Array data from the MOJAVE Program to get infoiTnation on the parsec-scale structure. 
Moreover, several Swift target of opportunity observations were triggered, adding important information on the source behaviour from optical/UV 
^ to hard X-rays. We paid particular attention to the calibration of the SwiftfUYCfT data to make it suitable to the blazars spectra. Simultaneous 
C observations from radio to y-rays allowed us to study in detail the correlation among the emission variability at different frequencies and to 
investigate the mechanisms at work during this high activity state of the source. 

In the period 9-30 March 2009, AGILE detected y-ray emission from PKS 1510-089 at a significance level of 21.5-(T with an average flux over 
• the entire period of (31 1 ± 21) x 10"* photons cm"^ s"' for photon energies above 100 MeV, and a peak level of (702 + 131) x 10"* photons cm"- 
'■ s"' on daily integration. The activity detected in y-rays occurred during a period of increasing activity from near-infrared to UV, as monitored 
by GASP-WEBT, REM and Swift/WOT. A flaring episode on 26-27 March 2009 was detected from near-IR to UV, suggesting that a single 
I mechanism is responsible for the flux enhancement observed at the end of March. By contrast, SwiftpiKF observations seem to show no clear 
I ' correlation of the X-ray fluxes with the optical and y-ray ones. However, the X-ray observations show a harder photon index (Fx- 1.3-1.6) with 
. . ' respect to most FSRQs and a hint of harder-when-brighter behaviour, indicating the possible presence of a second emission component at soft 
^ ' X-ray energies. Moreover, the broad band spectrum from radio-to-UV confirmed the evidence of thermal features in the optical/UV spectrum of 
k>( 1510-089 also during high y-ray state. On the other hand, during 25-26 March 2009 a flat spectrum in the optical/UV energy band was 

I^nJ . observed, suggesting an important contribution of the synchrotron emission in this part of the spectrum during the brightest y-ray flare, therefore 
a significant shift of the synchrotron peak. 
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1. Introduction 

According to the current paradigm (see Urry & Padovani 1995) 
Active Galactic Nuclei (AGNs) are galaxies whose emission is 
dominated by a bright central core, including a super massive 
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* The radio-to-optical data presented in this paper are stored in the 
GASP-WEBT archive; for questions regarding their availability, please 
contact the WEBT President Massimo Villata. 



black hole as central engine, surrounded by an accretion disk and 
by fast-moving clouds under the influence of the strong gravita- 
tional field, emitting Doppler broadened lines. Absorbing ma- 
terial in a flattened configuration, idealized as a dust torus, ob- 
scures the central parts so that for transverse line of sight only 
the more distant narrow-line emitting clouds are seen directly. In 
radio-loud objects we have the additional presence of a relativis- 
tic jet, roughly perpendicular to the disk. Within this scheme, 
blazars represent the fraction of AGNs with their jet at smaller 
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angles with respect to our line of sight, which causes relativistic 
aberration and emission ampUfication (Blandford & Rees 1978). 

Blazars are the most enigmatic subclass of AGNs, charac- 
terized by the emission of strong non-thermal radiation across 
the entire electromagnetic spectrum, from radio to very high y- 
ray energies. Because of their peculiar properties and amplified 
emission, blazars ofi'er the unique possibiUty to probe the cen- 
tral region of AGNs, and shed light on the mechanism respon- 
sible for the extraction of energy from the central black hole 
and for the acceleration and colhmation of relativistic electrons 
into jets. With the detection of several blazars in the y-rays by 
EGRET (Hartman et al. 1999) the study of this class of objects 
has made significant progress but, despite the big efibrts devoted 
to the investigation of the mechanisms responsible for the emis- 
sion in blazars, the definitive answer is still missing. The new 
generation of high-energy space missions, the Astrorivelatore 
Gamma a Immagini LEggew (AGILE) and Fermi Gamma-ray 
Space Telescope (GST) satellites, have to address some of the 
fundamental issues that were left unresolved by EGRET. 

PKS 15 10-089 is a nearby (z=0.361) blazar belonging to the 
class of the flat spectrum radio quasars (FSRQs) with radiative 
output dominated by the y-ray emission, while the synchrotron 
emission peaks around IR frequencies, below a pronunced UV 
bump, likely due to the thermal emission from the accretion disk 
(Malkan & Moore 1986; Pian & Treves 1993; D'Anomando et 
al. 2009a). 

Its radio emission exhibits very rapid, large amplitude vari- 
ations in both total and polarized flux (AUer, Aller, & Hughes 

1996) . Moreover, the radio jet shows superluminal motion with 
an apparent velocity higher than 20c, among the fastest of all 
blazars observed thus far (Jorstad et al. 2005), with the parsec 
and kiloparsec scale jets apparently misaligned by ~180 degrees 
(Homan et al. 2002, Wardle et al. 2005). 

PKS 1510-089 has been extensively observed in the X- 
rays during the last three decades, since the observations by the 
Einstein sateflite in 1980s (Canizares & White 1989). The X-ray 
spectrum as observed by ASCA in the 2-10 keV band was very 
hard with photon index of F ^ 1.3 (Singh, Shrader & George 

1997) , but steepened (T ^ 1.9) in tiie ROSAT bandpass (0.1- 
2.4 keV), suggesting the presence of a spectral break around 1-2 
keV, associated with the existence of a soft X-ray excess (Siebert 
et al. 1998). Subsequent observations by BeppoSAX (Tavecchio 
et al. 2000) and Chandra (Gambill et al. 2003) confirmed the 
presence of a soft X-ray excess below 1 keV. Evidence of a sim- 
ilar soft X-ray excess has been detected in other blazars such as 
3C 273, 3C 279, AO 0235-H164, and 3C 454.3. The origin of this 
excess is still an open issue, not only for blazars but for all AGNs 
(see e.g. D' Ammando et al. 2008, and the references therein). 

The interest in the peculiar X-ray spectrum of the source 
led up to a monitoring campaign organized during August 2006 
by Suzaku and Swift. The Suzaku observations confirmed the 
presence of a soft X-ray excess, suggesting that it could be 
a feature of the bulk Comptonization, whereas the SwiftfKKY 
observations revealed significant spectral evolution of the X- 
ray emission on timescales of a week: the X-ray spectrum be- 
comes harder as the source gets brighter (Kataoka et al. 2008). 
Unfortunately no y-ray satellites were operating at that time, 
therefore it was not possible to investigate the correlation be- 
tween X-ray and y-ray behaviour. However, y-ray emission from 
PKS 15 10-089 had aheady been detected in flie past by EGRET 
during low/intermediate states and was found to be only shghtly 
variable, with an integrated flux above 100 MeV varying be- 
tween (13 + 5) and (49 ± 18) x 10"^ photons cm"^ s"^ (Hartman 
etal. 1999). 



Instead, in the last three years PKS 1510-089 showed ex- 
treme variabihty over the whole electromagnetic spectrum and 
in particular very intense y-ray activity was detected, inter- 
mittently, by the AGILE and Fermi satellites. The monitoring 
of the source between September 2008 and June 2009 with 
the Large Area Telescope (LAT) onboard Fermi is summa- 
rized in Abdo et al. (2010a). The Ganmia-Ray Imaging Detector 
(GRID) onboard AGILE detected flaring episodes in August 
2007 (Pucefla et al. 2008) and March 2008 (D'Ammando et 
al. 2009a). Subsequently, an extraordinary y-ray activity was 
detected in March 2009, with several flaring episodes and a 
flux that reached ~ 700 X 10"^ photons cm"'^ s"^. Preliminary 
results were presented in D'Ammando et al. (2009b), Pucella 
et al. (2009), Verceflone et al. (2009). Recendy, Marscher et 
al. (2010) studying the flux behaviour at difl'erent frequencies, 
linear optical polarization, and parsec-scale structure of PKS 
1510-089 found a correspondence between the rotation of the 
optical polarization angle and the y-ray activity during the first 
six months of 2009. 

In this paper we discuss the results of the analysis of the multi- 
wavelength data of PKS 1510-089 coUected by GLAST-AGILE 
Support Program (GASP) of the Whole Earth Blazar Telescope 
(WEBT), Rapid Eye Mount (REM), Swift and AGILE during 
March 2009 and their implications for the emission mecha- 
nisms at work in this source. The paper is organized as fol- 
lows. Section 2 briefly introduces the multiwavelength cover- 
age on PKS 1510-089. In Sections 3 through 7 we present the 
analysis and results of AGILE, Swift, GASP- WEBT, REM, and 
Very Long Baseline Array (VLBA) data, respectively. In Section 
8 we discuss the contribution of the thermal emission in the 
low-energy part of the spectrum and the correlation between the 
emission at difi'erent energy bands. Finally, in Section 9 we draw 
our conclusions. 

Throughout this paper the quoted uncertainties are given at 
1-cr level, unless otherwise stated, and the photon indices are 
parameterized as A^(£') oc £ (ph cm s"' keV"' or MeV"' ) 
with r = or -I- 1 (or being the spectral index). 

2. The multiwavelength coverage 

The high y-ray activity observed by AGILE during March 2009 
triggered 15 Swift ToO observations, starting from 11 March 
2009, for a total of 46 ks. Moreover, the monitoring by GASP- 
WEB T in February-March 2009 provided important informa- 
tion from radio-to-optical band and, together with the data col- 
lected in near-IR and optical bands by REM during March 2009, 
allowed us to also obtain an excellent coverage in the low-energy 
part of the broad band spectrum of the source. Finally, high res- 
olution radio VLBI data at 15 GHz were obtained by monitoring 
the source within the Monitoring of Jets in Active galactic nuclei 
with VLBA Experiment (MOJAVE) project. 

A summary of the complete multiwavelength data set on 
PKS 1510-089 presented in this paper can be found in Table 
1. 



3. AGILE observations: data analysis and results 

The AGILE scientific instrument (Tavani et al. 2009) is very 
compact and combines four detectors that provide broad band 
coverage from hard X-rays to y-rays: a SiUcon Tracker (ST) op- 
timized for y-ray imaging in the 30 MeV - 30 GeV energy range 
(Prest et al. 2003; Barbiellini et al. 2001); a co-aligned coded- 
mask X-ray imager sensitive in the 18-60 keV energy range 
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Waveband 


Observatory 


Frequency/B and 


Radio 


SMA 


230 GHz 




Noto 


43 GHz 




Metsahovi 


37 GHz 




VLBA 


15, 43 GHz 




UMRAO 


8.0, 14.5 GHz 




Medicina 


5, 22 GHz 


Near-IR 


REM 


J, H, K 




Campo Imperatore 


J, H, K 




Roque (Liverpool) 


H 


Optical 


Abastumani 


R 




Calar Alto 


R 




Castelgrande 


B, V, R, I 




L'Ampolla 


R 




La Silla (MPG/ESO) 


B, V, R, I 




Lulin (SLT) 


R 




Roque (KVA and Liverpool) 


R 




San Pedro Martir 


R 




St. Petersburg 


B, V, R, I 




Valle d'Aosta 


R 




REM 


V, R, I 




Swift/UVCJT 


V, b, u 


UV 


Swift/UYOT 


uvwl, uvm2, uvw2 


X-ray 


SwiftPiKT 


0.3-10 keV 




Swift/BAT 


15-50 keV 


Gamma-ray 


AGILE GRID 


100 MeV - 30 GeV 



Table 1. Observatories contributing to the presented data set of 
PKS 1510-089 at different frequencies. 



(SuperAGILE; Feroci et al. 2007); a non-imaging CsI(Tl) Mini- 
Calorimiter sensitive in the 0.3-100 MeV energy range (MCAL; 
Labanti et al. 2009); and a segmented anticoincidence system 
(ACS; Perotti et al. 2006). The combination of ST, MCAL and 
ACS forms the GRID that assures the y-ray detection. 

The AGILE satellite observed PKS 1510-089 between 1 
March 2009 00:01 UT and 31 March 2009 11:41 UT (JD 
2454891.5-2454921.0), for a total of 325 hours of effective ex- 
posure time, during which the source moved from ~ 25° to ~ 55° 
off the AGILE pointing direction. 

AGILE-GRID data were analyzed, starting from the Level- 
1 data, using the AGILE Standard Analysis Pipeline and the 
AGILE Scientific Analysis Package. Counts, exposure, and 
Galactic background y-ray maps, the latter based on the dif- 
fuse emission model developed for AGILE (Giuliani et al. 2004), 
were generated with a bin size of 0.25° x 0.25° for photons with 
energies E > 100 MeV. We selected only the events flagged 
as confirmed y-ray events, and not collected during the South 
Atlantic Anomaly or whose reconstructed directions form an- 
gles with the satellite-Earth vector smaller than 80°, in order to 
reduce the y-ray Earth albedo contamination. 

PKS 1510-089 was detected over the period 1-30 March 
2009 at a significance level of 19.9-cr, with an average y-ray flux 
of /^E>iO()MeV = (162 + 12) X 10"*^ photons cm"^ s"'. Figure[T] 
shows the y-ray light curve of March 2009 with 1-day resolu- 
tion for photons with energy higher than 100 MeV. The average 
y-ray flux as well as the daily values were derived by the max- 
imum likelihood analysis, according to the procedure described 
in Mattox et al. (1993): first, the entire period was analyzed to 
determine the diffuse emission parameters, then the source flux 
density was estimated independently for each of the 1-day peri- 
ods with the diffuse parameters fixed at the value obtained in the 
first step. At the beginning of the observation period the source 
was not so active in y-rays as later AGILE-GRID does not detect 
PKS 1510-089 at a significance level higher than 3-cr between 



1 and 8 March 2009 and only upper limits at 95% confidence 
level are obtainecQ. Instead, in the period 9-30 March 2009 (JD 
2454899.5-2454921.5) different flaring episodes were detected 
and considering only this time interval the significance of detec- 
tion slightly increases to 21.5-cr, with an average flux of (31 1 + 
21) X 10"^ photons cm"^ s"'. 

In the period 9-30 March 2009 the AGILE 95% maximum 
likelihood contour level barycenterof the source is 1 = 351.278° 
b - 40.058°, with a distance between the y-ray position and the 
radio position (1 = 351.289°, b = 40.139°) of 0.08° and an overafl 
AGILE eiTor circle of radius r - 0. 15°. 

The y-ray light curve shows a complex structure with out- 
bursts having a duration of 4-5 days and approximately symmet- 
rical profile that could be an indication that the relevant timescale 
is the light crossing time of the emitting region. The peak level 
of activity with daily integration was Fe> lOOMeV = (702 ±131) 
X 10-** photons cm-2 s"' on 25 March 2009 (JD 2454916.0B 
This is the highest y-ray flux observed from this source on daily 
timescale and one of the highest fluxes detected from a blazar. 
We note that the increasing of the error on the flux estimation 
throughout the whole observation period is related to the fact 
that with the increase of the off-axis angle between the center 
of the field of view of the GRID and the position of the source, 
the possible background contamination and the GRID calibra- 
tion uncertainties become slightly higher Despite the different 
energy range and Tj^,, of the daily light curves, the AGILE re- 
sults are in agreement with the Fermi-LAT results presented in 
Abdo etal. (2010a). 

Unfortunately, the source was located substantially off-axis 
in the SuperAGILE field of view during the whole observation 
period, thus preventing us to extract a constraining upper limit 
in the 18-60 keV energy band with SuperAGILE. 

The y-ray spectrum during the first period of high activ- 
ity detected by AGILE (9-16 March 2009; JD 2454899.5- 
2454907.5) can be fitted with a power law of photon index F 
- 1.95 ± 0.15. The photon index is calculated with the weighted 
least squares method by considering four energy bins: 100-200 
MeV, 200-400 MeV, 400 MeV-1 GeV, and 1-3 GeV. A sim- 
ilar value of the photon index (F - 1.96 + 0.24) is estimated 
also in the period of the highest activity (23-27 March 2009; 
JD 2454913.5-2454918.5). Considering that during this second 
period the source was far off-axis with respect to the center of 
the GRID field of view, the energy bin 1-3 GeV is not consid- 
ered in order to reduce the possible background contamination 
for high off-axis angles at these energies. Comparing with the 
average photon index measured by EGRET (F - 2.47 + 0.21, 
Hartman et al. 1999) and that measured by Fermi-LAI over the 
first eleven months of observations (F - 2.41 + 0.01; Abdo et 
al. 2010b), the values observed by AGILE confirm a hardening 
of the y-ray spectrum during the very high activity detected in 
March 2009. A similar hard y-ray spectrum of PKS 1510-089 
was observed by AGILE during the flaring episode of March 
2008 (D'Ammando et al. 2009a). 

The coiTelation between the flux level and the spectral slope 
in the y-ray energy band was extensively studied by means of the 
analysis of all the EGRET data, but a firm result for the blazars 
is not found (Nandikotkur et al. 2007). A 'harder-when-brighter' 
behaviour in y rays of PKS 1510-089 was detected by Fermi- 
LAT over a long timescale (August 2008-June 2009) consid- 



' When the significance level of the detection obtained by the max- 
imum likelihood analysis is < 3-cr an upper limit at 95% confidence 
level is calculated (see Mattox et al. 1996). 

- Integrating between 0:00 UT and 23:59 UT for each day. 
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Fig. 1. AGILE y-ray light curve of PKS 1510-089 between 1 and 30 March 2009 (JD 2454891.5-2454921.5) for E > 100 MeV at 
1-day resolution. The downward arrows represent 2-cr upper limits. 



ering only flux levels of F (E > 200 MeV) > 2.4 x 10 pho- 
tons cnT^ s"' (Abdo et al. 2010a). A similar trend was observed 
in the EGRET era only for 3C 279 (Haitman et al. 2001) and 
marginally for PKS 0528+134 (Mukherjee et al. 1996), but this 
could be due to the fact that these were the only objects for which 
a long-term monitoring in y-rays was performed by EGRET. At 
present, with AGILE and Fermi we are able to follow a large 
number of blazars on very long timescales and investigate this 
behaviour on a larger sample of objects. In this context, consid- 
ering the long period of high activity recently observed in y-ray s, 
PKS 1510-089 is one of the best test cases, together with 3C 
454.3 (see e.g. Vercellone et al. 2010). 



4. Swift observations: data analysis and results 

The Swift satellite (Gehrels et al. 2004) performed 15 ToOs 
on PKS 1510-089 between 11 and 30 March 2009, triggered 
by the high y-ray activity of the source. The observations 
were performed with all the three onboard instruments: the X- 
ray Telescope (XRT; Buitows et al. 2005, 0.2-10.0 keV), the 
Ultraviolet Optical Telescope (UVOT; Roming et al. 2005, 170- 
600 nm), and the Burst Alert Telescope (BAT; Barthelmy et 
al. 2005, 15-150 keV). 



4.1. Swift/BAT data 

As part of its normal operations, BAT collects data over a wide 
area of the sky in its survey mode. The survey data in the 15-50 
keV band is used to produce sky images in which hard X-ray 
sources can be detected using the standard Swift analysis soft- 
ware, following the procedures described in Krimm et al. (2008, 
and the reference therein; see alscB). 

During March 2009, we detected two short flaring episodes 
from PKS 15 10-089. The first covered approximately 2 days be- 
ginning on 8 March 2009 (JD 2454899.0), with an average count 
rate of (0.006 + 0.002) counts s"' cm^^ (15-50 keV), which cor- 
responds to 28 mCrab and peaking on 9 March 2009 at 40 mCrab 
(Krimm et al. 2009). A second weaker episode occurred on 29 
March 2009 (JD 2454920.0), where the average count rate was 
(0.003 + 0.001) counts s"' cm"^, corresponding to 15 mCrab. 

As a comparison we report that during the period before 
the March 2008 y-ray flare detected by AGILE (D'Ammando 



http://swift.gsfc.nasa.gov/docs/swift/results/transients 



et al. 2009a), Swift/BAil recorded flaring episodes from PKS 
1510-089 on 16 Febmary 2008 (JD 2454513.0), 21 February 
2008 (JD 2454518.0) and 24 February 2008 (JD 2454521.0), 
with count rates of approximately (0.005 ± 0.002) counts s"' 
cm"^ on each of those observations. 



4.2. Swift/XRT data 

The XRT data were processed with standard procedures 
(xrtpipeline vO.12.4), filtering, and screening criteria by us- 
ing the Heasof t package (v.6.8). The source count rate was low 
during the whole campaign (count rate < 0.5 counts s '), so 
we only considered photon counting (PC) data and further se- 
lected XRT event grades 0-12. Pile-up correction was not re- 
quired. Source events were extracted from a circular region with 
a radius between 15 and 20 pixels (1 pixel ~ 2.36"), depend- 
ing on the source count rate, while background events were ex- 
tracted from a circular region with radius 40 pixels away from 
background sources. Ancillary response files were generated 
with xrtmkarf, and account for different extraction regions, vi- 
gnetting and PSF corrections. We used the spectral redistribu- 
tion matrices vOll in the Calibration Database maintained by 
HEASARC. Considering the very low number of photons col- 
lected on 26 and 27 March we summed the data of these two 
consecutive days in order to have enough statistics to obtain a 
good spectral fit. 

All spectra were rebinned with a minimum of 20 counts per 
energy bin to allow fitting within XSPEC (vl2.5.1; Amaud 
1996). We fit the individual spectra with a simple absorbed 
power law, with a neutral hydrogen column fixed to its Galactic 
value (6.89 x 10^" cm^^; Kalberla et al. 2005). The fit results 
are reported in Table |2] The comparison between our results and 
those presented by Abdo et al. (2010a) reveals some differences, 
which we ascribe to their using a lower number of counts per 
bin than generally adopted to allow fitting (> 20 counts/bin). 
Indeed, these differences increase when the number of counts 
decreases. 

During the 15 ToOs performed in March 2009, SwifiPiKT 
observed the source with a 0.3-10.0 keV flux in the range (0.8 - 
1.1) X 10"" erg cm"^ s"', a somewhat fainter state with respect 
to the high flux level observed in August 2006 (FIuxq 5_iokeV 
- (1.1 - 1.8) X 10"" erg cm"^ s"'). Instead, a hint of spectral 
trend seems to be present in the X-ray data. Figure |2] shows the 
XRT photon indexes as a function of the X-ray fluxes in the 0.3- 
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Table 2. Log and fitting results of Swift/XRT obsei-vations of PKS 1510-089 during March 2009. Power law model with A^h fixed 
to Galactic absorption (6.89 x 10^*^ cm"-') is used. " Unabsorbed flux. 



Start Time 


Stop Time 


Exp. Time 


Photon Index 


Flux 0.3- 


10.0 keV 


Xli (d-o.f ) 


(yyyy-mm-dd hh:mm:ss) 


(yyyy-mm-dd hh:mm:ss) 


(sec) 




F 


(X 10"'^ er 


g cm"^ s"' 


) 




9000-01-11 91-9'S-'Sfi 


4801 
^oy i 


1 


_i_ n 1 n 


9.28 ±0.77 


Q1Q ns^ 


9000-0^-19 1^-1 ^-4.4. 


9000-01-19 90-90-'S0 


4842 


i.4D 


± u.uy 


10.86 


±0.84 


U.OU / \JO } 


900Q-0^-17 04-9Q-17 


9000-01-17 11 •07-^6 


48(S0 


1 ^ 1 
i.ji 




8.67- 


tO.77 


1 041 ni^ 


900Q-0^-ix oo-n9-in 
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10 keV: a possible hardening of the spectrum with the increase 
of the flux is observed, in spite of the large errors. A similar 
harder- when-brighter behaviour in X-rays was observed for PKS 
1510-089 in August 2006 (Kataoka et al. 2008), March 2008 
(D'Ammando et al. 2009a), and over the period January-June 
2009 (Abdo et al. 2010a). 



The significant spectral evolution observed by XRT, with the 
photon index changing from 1.3 to 1.6, is not usual for FSRQs, 
for which only little X-ray variability is detected both on short 
and long timescales. Moreover, the average X-ray photon index 
observed over March 2009 for this object (< F > = 1.48 ± 0.03) 
is lower than that of the radio-loud quasars (< F > = 1.66 + 
0.07; Lawson et al. 1992, Cappi et al. 1997) and is more similar 
to that of the high-redshift (z > 2) quasars (e.g. Page et al. 2005). 
Since usually in FSRQs the X-ray energy range samples the low- 
energy tail of the external Compton (EC) component, as already 
discussed in Kataoka et al. (2008) and Abdo et al. (2010a), such 
hard spectral indexes in X-rays should imply a very flat energy 
distribution, challenging the standard shock models of particle 
acceleration. In a standard shock model these hard X-ray spec- 
tra produce relativistic electrons with an energy spectrum harder 
than the canonical power law distribution N(y) oc y^^, or alter- 
natively, as discussed in Sikora et al. (2002), another mecha- 
nism (e.g. instabilities driven by shock-reflected ions, Hoshino 
et al. 2002; or magnetic reconnection, Romanova & Lovelace 
1992) should energize the low energies electrons that typically 
produce the X-ray emission in the EC model. 



4.3. Swift/UVOT data 



UVOT observed PKS 1510-089 in all its optical (v, b, and u) 
and UV (uvwl, uvml, and uvwl) photometric bands. Data were 
reduced with the uvotmaghist task of the HEASOFT package. 
Source counts were extracted from a circular region of 5 arcsec 
radius, centred on the source, while the background was esti- 
mated from a surrounding annulus with 8 and 18 arcsec radii. In 
Fig.[3]UVOT magnitudes are displayed with blue circles. 
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Fig. 2. SwiftpiKT photon index as a function of the 0.3-10 keV 
flux. 



5. GASP observations in radio, near-IR and optical 
bands: data analysis and results 

The GLAST-AGILE Support Program (Villata et al. 2008, 
2009a) is a project bom from the Whole Earth Blazar Telescop^ 
in 2007. It is aimed at providing long-term monitoring in the op- 
tical (/?-band), near-IR, and radio bands of 28 y-ray-loud blazars 
during the lifetime of the two y-ray satellites. 

The /;-band GASP observations of PKS 1510-089 in the 
period considered in this paper were performed by the fol- 
lowing observatories: Abastumani, Calar Altc0, Castelgrande, 
L'AmpoUa, La Sifla (MPG/ESO), Lulin (SET), Roque de 
los Muchachos (KVA and Liverpool), San Pedro Martir, St. 
Petersburg, and Valle d'Aosta. These data were calibrated with 
respect to stars 2-6 by Raiteri et al. (1998). The GASP obser- 
vation in /?-band showed that after a low intensity period in 
February 2009, with the source observed constantly at ~ 16.5 



* http : //www . oato . inaf . it/blazars/webt J. 

' Calar Alto data were acquired as part of the MAPCAT project 
|http://www.iaa.es/~ia gudo/research/ MAPCAT( 
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Fig. 3. Light curves collected in near-IR, optical and 
UV bands, between 5 Mai'ch and 2 April 2009 
(JD 2454895.5-2454924.5). Blue circles represent the 
UVOT data in v, b, u, uvwl, uvml, uvw2 filters. Red di- 
amonds represent GASP data in R, J, H, and K bands. 
Magenta triangles represent WEBT data in /, V, and B 
bands. REM data in V, R, I, J, H, and K bands are rep- 
resented with black squares. Yellow regions in i?-band 
light curve indicate the peaks of y-ray activity observed 
by AGILE. 



mag, the optical activity of the source greatly increased in March 
2009 (Villata et al. 2009b) with an intense flare on 27 March (JD 
~ 2454917.6), after a brightening by almost 1 mag in 3 days, 
reaching R = 14.35 + 0.03 mag (Larionov et al. 2009; see Fig. 3 
and Fig. 4). Near-IR data in J, H, and K bands are from Campo 
Imperatore and Roque de los Muchachos (Liverpool), whereas 
WEBT data in B, V, and / bands were taken at Castelgrande, 
La SiUa and St. Petersburg. Data collected by GASP and WEBT 
observatories during March 2009 are reported in Fig. |3]with red 
diamonds and magenta triangles, respectively. 

Radio fluxes were measured at: Submillimeter AiTay0(SMA, 
230 GHz), Medicina (22 and 5 GHz), Metsiihovi (37 GHz), Noto 



^ These data were obtained as part of the normal monitoring program 
initiated by the SMA (see Gurwell et al. 2007). 



(43 GHz), and UMRAO (14.5 and 8.0 GHz). The radio data are 
shown in Fig. ID together with the optical (/?-band) and near- 
IR (//-band) Hght curve of PKS 1510-089 in Febmary-March 
2009. 

We point out that in those cases where diff'erent datasets were 
present in the same band, we performed a careful data analysis 
to determine possible offsets and corrected for them to get ho- 
mogeneous light curves, as usually done for GASP- WEBT data 
(see e.g. Viflata et al. 2002, Raiteri et al. 2005). 

6. REM observations: data analysis and results 

During March 2009 photometric near-IR and optical observa- 
tions were carried out with REM (Zerbi et al. 2004), a robotic 
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Fig. 4. ;?-band light curve of PKS 15 10-089 obtained by GASP 
and REM during February-March 2009, together with the H 
band data by GASP (top panel), compared with the radio flux 
densities at different frequencies (bottom panel). 



stationary, and three emitted knots. From linear regression fits 
we obtain highly superluminal apparent separation velocities 
between 16c and 20c, consistent with the separation speed 
derived for other knots (e.g. Homan et al. 2002), and at different 
frequency (Marscher et al. 2010), indicating an intrinsic separa- 
tion velocity /Sinn- > 0.9989 and an orientation with the line of 
sight < 5°. 

From the pc-scale resolution images we could separate the 
core flux density from the jet emission, in order to study 
the lightcurve and polarization properties of each of them. 
Between 2007 and 2010, three main episodes of enhanced 
radio luminosity have been detected, successively followed by 
marked flux density decreases. By contrast, no obvious trend of 
the polarization percentage has been found. 
After March 2009, the flux density of the core region, that was 
in a minimum, started to increase, while the flux density of the 
jet component did not show any significant luminosity changes. 
A remarkable aspect observed just after the March 2009 flare 
is the abrupt change of about 75° in the polarization angle of 
the core, whereas the fractional polarization remains constant. 
More details on the long-term radio monitoring of this source 
are presented in Orienti et al. (2010). 



telescope located at the ESO Cerro La Silla observatory (Chile). 
The REM telescope has a Ritchey-Chretien configuration with 
a 60 cm f/2.2 primary and an overall f/8 focal ratio in a fast 
moving alt-azimuth mount providing two stable Nasmyth focal 
stations. At one of the two foci, the telescope simultaneously 
feeds, by means of a dichroic, two cameras: REMIR for the near- 
IR (Conconi et al. 2004) and ROSS (Tosti et al. 2004) for the 
optical. Both cameras have a field of view of 10 x 10 arcmin 
and imaging capabilities with the usual near-IR (z, J, H, and 
K) and Johnson-Cousins VRI filters. The REM software system 
(Covino et al. 2004) is able to manage complex observational 
strategies in a fully autonomous way. All raw near-IR/optical 
frames obtained with the REM telescope were then corrected for 
dark, bias, and flat field following standard recipes. Instrumental 
magnitudes were obtained via aperture photometry and absolute 
calibration has been performed by means of secondary standard 
stars in the field (see Raiteri et al. 1998; Smith and Balonek 
1998). The data presented here were obtained by GO program 
for AOT19 (PI: F D'Ammando). REM data in V, R, I, J, H, and 
K bands are reported in Fig. |3]with black squares. 

7. VLBA observations: data analysis and results 

When observed with the high spatial resolution provided by 
the Very Long Baseline Interferometer (VLBI) technique, 
PKS 1510-089 shows a radio structure dominated by the core 
region from which the pc-scale jet emerges forming an angle 
of -28°, i.e. in the north-west direction. Monitoring campaigns 
pointed out the ejection of new superluminal components, and 
followed their evolution throughout several observing epochs 
(see Marscher et al. 2010). 

As part of the MOJAVE program (Lister et al. 2009), PKS 
1510-089 is frequently observed with the VLBA at 15 GHz 
(see Fig. |5]l. To study variations in the source structure, we 
retrived the Mv-datasets and we imported them into NRAO AIPS 
package. Final images both in total intensity and polarization 
have been produced after a few phase-only self-calibration 
iterations (for more details on calibration and image analysis 
see Orienti et al. 2010). From the analysis of the multi-epoch 
observations obtained between 2007 and 2010 it was possible 
to follow the separation between the core region, considered 
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Fig. 5. VLBA image at 15 GHz of PKS 1510-089 on 25 March 
2009. On the image we provide the restoring beam, plotted in 
the bottom left corner, the peak flux density in mJy beam and 
the first contour (f.c.) intensity in mJy beam"', which is 3 times 
the off-source noise level. Contour levels increase by a factor of 
2. 
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8. Discussion 

8.1. Thermal signatures in tlie radio-to-UV spectrum 

Evidence of thermal emission in blazars has been detected in a 
limited number of objects until now, and usually observed during 
low activity states (e.g. 3C 279, Pian et al. 1999; 3C 273, Grandi 
& Palumbo 2004; 3C 454.3, Raiteri et al. 2007; AO 0235+164, 
Raiteri et al. 2008; NRAO 150, Acosta-Pulido et al. 2010). This 
is because, contrary to the other AGNs, in blazars the jet emis- 
sion at small angles with respect to the line of sight of the ob- 
server is strongly amplified due to the relativistic beaming ef- 
fect, overwhelming all the thermal contributions. However, these 
thermal features are less prominent in blazars but not absent, at 
least for FSRQs, which present a radiatively efficient accretion 
disk and broad emission lines. 

The knowledge of these thermal features could have im- 
portant consequences not only for the low-energy part of the 
spectrum, in which the thermal emission could be directly ob- 
served, but also for the high-energy part of the spectrum, because 
the photons produced by the accretion disk, directly or through 
the reprocessing by the broad line region (BLR) or the dusty 
torus, can become the main source of seed photons for the EC 
mechanism that usually dominates the y-ray emission of FSRQs 
(Ghisellini et al. 1998). 

In this context PKS 1510-089 seems to be peculiar. Since 
the synchrotron peak of this source is usually observed in mid- 
infrared band (see Bach et al. 2007; Nieppola et al. 2008), the 
optical/UV part of the spectrum is not dominated by the non- 
thermal jet emission, allowing us to observe directly the ther- 
mal manifestation of both the accretion disk (the so-called "big 
blue bump", e.g. Laor 1990) and of the BLR (the so-called "lit- 
tle blue bump", e.g. Wills et al. 1985) also during high activity 
states, like in mid-March 2008 (see D'Ammando et al. 2009a). 
Moreover, an excess at far-IR wavelengths was observed in PKS 
1510-089 by IRAS (Tanner et al. 1996), likely due to dust ra- 
diation from the torus. Recently, Spitzer observations of PKS 
1510-089 produced only an upper limit of 2 x 10"^^ erg s"' to 
the thermal emission from dust torus (Malmrose et al. 201 1). 

In order to distinguish between non-thermal and ther- 
mal emission contributions, we built spectral energy distribu- 
tions (SEDs) from radio to UV with GASP-WEBT, REM and 
SwiftfXJVOT data. Conversion of magnitudes into de-reddened 
flux densities was obtained by adopting the Galactic absorption 
value Ab = 0.416 from Schlegel et al. (1998), the extinction laws 
by Cardelli et al. (1989) and the magnitude-flux calibrations by 
Bessefl et al. (1998). 

As for the 5w//f/UVOT data, we noticed that PKS 1510-089 
has a b-v ~ 0.3 that is out of the validity range indicated by 
Poole et al. (2008) for their flux calibrations in the UV bands. 
Hence, following Raiteri et al. (2010, 201 1) we calculated efi'ec- 
tive wavelengths /leff, count rate to flux density conversion fac- 
tors CF, and amount of Galactic extinction Aa for each UVOT 
band, by folding the quantities of interest with the source spec- 
trum and effective areas of UVOT filters. The results are shown 
in Table 3. 

UVOT count rates were thus converted into flux densities 
and then corrected for Galactic extinction adopting the CF and 
A A reported in Table 3. In Figure 6 we display three SEDs of 
March 2009 as weU as the March 2008 SED already published 
by D'Ammando et al. (2009a). The comparison between the 
March 2008 SED obtained with standard calibrations (Figure 3 
of D'Ammando et al. 2009a) and the same SED obtained by 
applying the procedure described above reveals that the dip in 
the uvwl band remains, though it is now less pronounced. This 



Table 3. Results of flie UVOT cahbration procedure for PKS 
1510-089. 



Filter 




CF 


Aa 




A 


10""'ergcm"2 s"' A"' 


mag 


V 


5422 


2.61 


0.33 


b 


4346 


1.47 


0.43 


u 


3466 


1.65 


0.52 


uvwl 


2633 


4.21 


0.73 


uvml 


2251 


8.45 


0.90 


uvwl 


2059 


6.31 


0.86 



means that it is not due to calibration problems, but likely re- 
flects a true intrinsic feature. Indeed, we expect the little blue 
bump due to BLR emission (Mgll, Fell, and Balmer lines) to 
peak in the u-b frequency range. Near-IR to UV spectra of PKS 
1510-089 in March 2009 have also been presented by Abdo et 
al. (2010a), but for a different choice of epochs, which prevents 
a detailed comparison with our results. However, as expected, 
we notice a deeper wl dip than that we obtained with our more 
accurate calibrations. 

In other FSRQs, such as 3C 454.3, the presence of the little 
and big blue bumps was detected only during low activity state 
of the source (see Raiteri et al. 2008, Vercellone et al. 2010). 
Instead, the SED of PKS 1510-089 collected on 18 March 2009 
confirmed the evidence of thermal signatures in the optical/UV 
spectrum of the source also during high y-ray states, with a con- 
tribution in the optical part likely due to the little blue bump 
and a significant rise of the spectrum at UV due to the accretion 
disk emission, as already observed on 20-22 March 2008. On 
the other hand, the broad band spectrum of PKS 15 10-089 from 
radio-to-UV during 25-26 March 2009 shows a flat spectrum 
in the optical/UV energy band, suggesting an important contri- 
bution of the synchrotron emission in this part of the spectrum 
during the brightest y-ray flaring episode and therefore a sig- 
nificant shift of the synchrotron peak, usually observed in this 
source in the infrared band. We derived an estimate of the fre- 
quency of the synchrotron peak in the three SEDs applying a 



T 




log V [Hz] 

Fig. 6. SEDs of the low-energy part of the spectrum constructed 
with data collected by GASP-WEBT (empty circle), KVA (plus 
sign), Abastumani (cross), St. Petersburg (asteriks), Campo 
Imperatore (diamonds), REM (squares), and Swift/UYOT (filled 
circles) during March 2008 and March 2009. A cubic polyno- 
mial fit was applied to the SEDs for locating the position of the 
synchrotron peak. 
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cubic polynomial fit (see e.g. Kubo et al. 1998) to the radio and 
IR data, which are likely due to pure synchrotron, as shown in 
Fig.|6] The synchrotron peak shifted from v = 1.5 x 10'^ Hz to v 
= 6.5 X 10'^ Hz between 18 and 26 March 2009. Similarly to the 
harder- when-brighterbehaviour observed in X-rays, this is a typ- 
ical behaviour of high-frequency-peaked BL Lacs (HBLs) and 
not so commonly observed in FSRQs such as PKS 1510-089, 
even if this could be partially due to the fact that it is more diffi- 
cult to obtain a long term monitoring of the synchrotron peak of 
FSRQs in the IR band with respect to the optical/UV and X-ray 
band. 

Abdo et al. (2010a) analyzing the UVOT data of PKS 
1510-089 with standard calibrations (and thus with possible bi- 
ases in the results, see above) found an anticorrelation between 
UV flux and UV hardness ratio in the period January-June 2009, 
suggesting a different level of contamination in UV by the high- 
energy branch of synchrotron emission depending on the activity 
level. This is in agreement with our results, even if the impor- 
tance of the increase of the synchrotron emission observed at the 
end of March 2009 seems to be much higher with respect to the 
general trend. 

The UVOT data reported in the SED of 25-26 March 2009 
are collected at JD = 2454916.46, instead the BVRI data are 
collected by St. Petersburg at JD - 2454916.51-53, thus the sep- 
aration in time between UVOT and BVRI data is about 1 hour 
This rules out the possible bias related to the optical/UV vari- 
ability of the object, confirming the flat optical/UV spectrum. 
The NIR data for the same SED are collected by REM at JD - 
2454916.66, about 3.5 hrs after the optical data and 4.5 hrs after 
the UVOT data. In this case it is not possible to completely rule 
out a mismatch due to the possible rapid optical/UV variabil- 
ity, but by comparison with the NIR data collected during 20-22 
March 2008 and 18 March 2009, it is evident the change of the 
NIR spectrum, in agreement with what we observed also in op- 
tical/UV. Moreover the comparison between the optical/UV data 
coUected by UVOT on 25 March 2009, during the first UVOT 
exposure (JD ~ 2454916.1, orange squares in Fig.|6) and those 
taken about 9-10 hours after (JD ~ 2454916.5, green triangles), 
reveals a noticeable spectral and flux variation. Considering that 
the accretion disk is slowly variable on such short timescales, 
this is another proof of the fact that the rapid increase of the op- 
tical emission started on 25 March and peaking on 27 March is 
mainly due to the synchrotron mechanism. This is also in full 
agreement with the simultaneous rapid increase of the degree of 
optical polarization shown in Marscher et al. 2010 (see in partic- 
ular their Fig. 4). Thermal emission is unpolarized as it reflects 
the random walk of atoms and ions within the emitting region, 
therefore such increase of the degree of polarization is a clear 
signature of a rapid and strong increase of the contribution in 
the optical band of a non-thermal mechanism such as the syn- 
chrotron emission. 

8.2. Light curve behaviour and correiations 

By comparing the source behaviour observed in y-rays by 
AGILE with the ones observed from optical-to-UV by GASP- 
WEBT, REM and Swift/UYOT during March 2009 we noted 
that while the y-ray light curve shows evidence of different out- 
bursts of increasing entity, the optical and UV light curves seem 
to show a gradual increase of the flux in time with a rapid flux 
enhancement after JD ~ 2454915.0 and a single major outburst 
occun-ed between JD 2454917.25 and JD 2454917.60. A varia- 
tion of 0.48 mag in 9 days (JD 2454906-2454915) was detected 
in the /?-band light curve, followed by a more rapid variation of 



0.95 mag in ~ 2 days. The uvwl light curve shows a similar be- 
haviour, with a variation of 0.43 mag in 7 days and 0.78 mag 
in the following 2 days. The larger increase observed in /?-band 
with respect to uvwl band is in agreement with a larger con- 
tribution of the synchrotron emission in the optical than in the 
UV band. A linear fit applied to the /?-band and Mvwl-band light 
curves in Fig. 3 shows a change of slope after JD = 2454915. 
The optical/UV peak seems to be delayed with respect to the 
brightest y-ray peak by 1-2 days, although we cannot exclude a 
very rapid optical/UV flare occurred during the gap in the optical 
and UV light curves before JD 2454917.25. 

Another possibility is that the delay between optical and y- 
ray emission is due to our choice of the To used for building 
the y-ray daily light curve. To investigate the influence of the Tq 
on the determination of the y-ray peak we constructed three dif- 
ferent light curves, considering an integration interval of 1 day 
and shifting the Tq of AT] - 8 hours and AT2 =16 hours with 
respect to our initial choice (Tq - 0:00 UT). In Fig. |7]the op- 
tical (R-hand), UV (wl filter) and the three y-ray light curves 
built with the different To in the period 21-31 March 2009 (la- 
belled as I, II and III) are shown. For each y-ray light curve we 
assign a weight "1" to the three fractional 8-hours bins that con- 
stitute the total 1-day bin with the highest flux and a weight "0" 
to the other fractional bins. Combining the information derived 
from the three light curves we can build a probability function 
and estimate the fractional bin that corresponds with the highest 
probability to the emission peak. In this way we estimated the 
y-ray emission peak at JD - 2454916.33 + 1.00, leading to a 
reduced delay between y-ray and optical/UV emission, in accor- 
dance also with the delay between the FermZ-LAT and optical 
light curve peaks showed in Abdo et al. (2010a) for the same pe- 
riod. Moreover, we note that a recent investigation of the Fermi- 
LAT data of PKS 1510-089 coflected in the first half of 2009 
by Tavecchio et al. (2010) showed significant y-ray variability 
on timescales of 6 hours and a flux peak on JD = 2454917.25, 
compatible with our result. 

However, a significant increase of the optical/UV emission 
is detected only for the brightest y-ray flare, whereas the opti- 
cal/UV activity before JD ~ 2454915 was increasing but only 
slowly, according to a dominant contribution of the thermal 
emission, over which a subsequent rapid variability of the syn- 
chrotron emission is superimposed (see Sec. 8.1). A progressive 
increase of the source activity, similar to the behaviour of the op- 
tical light curve collected by GASP-WEBT, was observed also 
in the near-IR and optical bands by the REM telescope (Fig.O, 
suggesting that the synchrotron is the dominant mechanism re- 
sponsible for the flux enhancement observed from near-IR to 
UV at the end of March 2009. In this context, the lag of 13+1 
days between the y-ray and the optical peaks reported in Abdo 
et al. (2010a) could suggest a different origin of the dominant 
component in optical band during the January and April 2009 
flares with respect to the flare occurred at the end of March 2009. 
This is another clue in favour of a more complex correlation be- 
tween the optical and y-ray emission in the FSRQs, in which 
the optical emission could be due to the contribution of both the 
synchrotron emission and the thermal accretion disk emission. 
Despite the different dominant contribution in the optical part of 
the spectrum, a similar photon index was observed by AGILE in 
the y-ray spectrum during the first activity period (9-16 March) 
and the highest activity period (23-27 March), suggesting that 
the same radiation mechanism is responsible for the high energy 
emission of the source over the whole observing period. 
Moreover, it is worth noting that the variability amplitude ob- 
served at the end of March 2009 is similar in optical (a factor of 
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Fig. 7. Comparison of the /?-band light curve of PKS 15 10-089 
(panel (a); GASP-WEBT: red squares, REM: black dots) with 
the uvwl band data by Swift/WOT (panel (b)), the X-ray data by 
SwiftpCRT (panel (c)), and the AGILE y-ray light curves in the 
three different configurations (panel (d), (e), (f), see Sec. 8.2 for 
details) during 21-31 March 2009 (JD 2454913.5-2454922.0). 



~3) and y rays (a factor of ~3-4). If the y rays were produced 
by inverse Compton scattering of photons reprocessed by the 
BLR, as suggested by Abdo et al. (2010a), the y-ray emission 
would be proportional to the number of the emitting electrons, 
like the synchrotron flux. Therefore, the nearly linear correlation 
observed between optical and y rays may be due to the variation 
of the electron number in the jet. In this context, the lag of 1 day 
of the optical peak with respect to the y-ray one could be only 
an upper limit due to the lack of data before the observed optical 
peak. 

Compared to the optical and y-ray activity, the X-ray flux is 
not very variable and does not seem to be correlated with the 
high y-ray activity (see Fig. 7, panel c). The lack of evident X- 
ray/y-ray correlation is not so surprising. In fact, Marscher et 
al. (2010) akeady reported that the long-term X-ray light curve 
of PKS 1510-089 collected during the 2006-2009 period by 
Rossi X-ray Timing Explorer is significantly correlated with the 
14.5 GHz one rather than with the y-ray emission. This evidence 
could be justified by the fact that the X-ray photons are likely 
to originate mostly in the low-energy part of the EC emission. 
Therefore the X-ray spectral index seems to reflect the change in 
the low-energy tail of the electron distribution, whereas a change 
in the high-energy part of the electron distribution should be the 
driver of the variation of the y-ray flux, as already pointed out by 



Abdo et al. (2010a). Moreover, as already discussed in Kataoka 
et al. (2008) and D'Ammando et al. (2009a), the spectral evolu- 
tion of PKS 1510-089 could be justified by the contamination 
at low energies of a second emission component with respect 
to the EC emission. In this case the X-ray spectral shape of the 
EC component should remain almost constant, regardless of the 
activity level, but its flux should vary. Therefore the amount of 
contamination of the second component should change, being 
more important when the source is fainter and almost hidden by 
the EC component when the source is brighter Possible origins 
proposed for this second component are the soft X-ray excess, 
the bulk Comptonization or a more significant contribution of 
the SSC component in any activity states. A possible hint of in- 
crease of activity in X-rays is detected by XRT on 28 March and 
this could be an indication of an enhancement of the contribution 
of the synchrotron self Compton (SSC) emission, strictly related 
to the optical/UV outburst observed on 27 March. The strength 
of the SSC with respect to EC emission depends on the ratio be- 
tween the synchrotron and the external radiation energy density, 
as measured in the comoving frame, U'^^JU'^^^. Taking into ac- 
count that, under the assumption that the dissipation region is 
within the BLR, U'^y^ depends on the injected power, the size 
of the emission and the magnetic field, a larger magnetic field 
can lead to a larger SSC contribution together with an increase 
and a shift of the synchrotron peak. By contrast, an increase of 
the SSC emission should lead to a not negligible contribution of 
this component, which usually shows a photon index of F ~ 2 
in X-rays for the FSRQs, and therefore a softening of the total 
X-ray spectrum. Instead a hard X-ray spectrum was observed by 
XRT (F.v = 1.31 + 0.14 on 28 March), confirming that also dur- 
ing a contemporaneous synchrotron flare the SSC contribution 
is not significant in the X-ray spectrum of PKS 1510-089, as 
already observed by Suzaku during a low state of the source in 
2009 (Abdo et al. 2010c). 

An alternative scenario is that the electrons radiating at 14.5 
GHz produce X-rays via IC scattering of the IR photons from 
a dust torus. This interpretation favours the "far dissipation" 
scenario recently proposed again for the emission of FSRQs 
by Sikora et. (2008) and Marscher et al. (2010), in which the 
high energy peak is mainly due to the EC scattering of the IR 
seed photons produced by the torus. In particular, Marscher et 
al. (2010) identified the dissipation zone with a moving superlu- 
minal knot before and after passing the core, which is interpreted 
as the first recollimation shock visible at millimeter wavelengths 
in the jet. 

It is crucial to discriminate between these two scenarios in 
order to understand where the plasma blob dissipates in the jet. 
In fact, the energy density of the different sources of seed pho- 
tons varies with the distance from the central black hole, there- 
fore the dominant population of target photons depends on the 
site of the dissipation region that produces the flaring episode. 
The determination of the location of the emitting region in the 
jet of blazars is an open issue, in particular for FSRQs. However, 
the rapid variability observed in the y-ray light curve of PKS 
1510-089 collected by AGILE in March 2009 (even more re- 
markable in the Ferm/-LAT light curves on timescales of hours, 
Abdo et al. (2010a), Tavecchio et al. 2010) puts severe con- 
straints on the size of the emitting region, leading to a very 
compact region difficult to be explained if the y-ray emission 
is produced far away from the central black hole. This seems 
to disfavour the "far dissipation" scenario. On the other hand, 
the recent detection of this source at TeV energies by H.E.S.S. 
(Wagner and Behera 2010) suggests possible complex scenarios. 
We can envisage that more than one emission region is simulta- 
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neously active and located at different distance from the central 
black hole in a spine-sheath (Tavecchio and Ghisellini 2008) or 
jet-in-jet structure (Giannios et al. 2009). 

It is worth noting that the outburst detected by ^vvZ/f/BAT 
on 8-9 March 2009 occurred just at the beginning of the y-ray 
activity observed by AGILE. On the other hand, no significant 
activity is detected by BAT simultaneously with the y-ray flares 
detected by AGILE. 

Finally, as shown in Figure H] no clear sign of strong activity 
in radio bands is observed during March 2009, although the radio 
coverage is sparse. However, the SMA data show an increase 
of about a factor of two between February and March 2009 at 
millimeter wavelengths, similarly to the optical behaviour even 
if with a minor variation amplitude. Moreover, the radio data 
acquired at Metsahovi at 37 GHz on a longer timescale show 
that the increase of activity marginally observed also in Fig. HI 
which starts at the beginning of March 2009 (JD ~ 2454910), 
has continued gradually in April-May 2009, reaching a peak flux 
density of 4.02 Jy on 15 May 2009 (JD 2454969.4; see Abdo et 
al. 2010a). The increase of flux density observed at 37 and 230 
GHz confirms that the mechanism producing the y-ray flaring 
events also interested the mm/cm emitting region with difference 
in time likely related to opacity effects, as already observed in 
March-April 2008 during the period soon after the previous y- 
ray flare of PKS 15 10-089 detected by AGILE (D'Ammando et 
al. 2009a). 

In addition, the VLBA data show a minimum of the flux den- 
sity of the core in March 2009 and since the rotation of the po- 
larization angle of the core was very mild in March relative to 
next 50 days, it is unlikely that the March y-ray flaring episodes 
are related to the the blob associated with later flares as in the 
Marscher et al. (2010) model. 

9. Summary and concluding remarks 

We reported on the detection by AGILE-GRID of y-ray activity 
from the FSRQ PKS 15 10-089 during March 2009. AGILE data 
showed an amazing y-ray activity, with several flaring episodes 
and a flux that reached on daily integration ~ 700 x 10 photons 
cm"^ s the highest flux detected from this source until now 
and one of the highest detected from a blazar. No significant 
spectral changes are observed between the different y-ray flares 
in March 2009, indicating that the same mechanism is dominant 
for the high energy emission of the source over all the period. 

This y-ray activity triggered simultaneous multiwavelength 
observations, which provided us a wide dataset for studying the 
correlation between the emission properties from radio to y rays. 
An increasing activity in near-IR and optical was observed, with 
a strong flux enhancement after 24 March and a flaring episode 
on 26-27 March 2009 observed by GASP-WEBT and REM, al- 
most simultaneous to the brightest y-ray flare. By contrast, no 
clear signs of simultaneous high activity were detected from 
near-IR to UV during the other y-ray flares observed by AGILE. 

The broad band spectrum from radio-to-UV collected on 18 
March confirmed the evidence of thermal signatures, the lit- 
tle and big blue bumps, in the optical/UV spectrum of PKS 
1510-089 also during high y-ray states, and not only during 
low activity states as observed in other blazars. The radio-to- 
UV spectrum on 25-26 March showed a flat spectrum in the 
optical/UV energy band, suggesting an important contribution 
of the synchrotron emission in this part of the spectrum during 
the brightest y-ray flare and therefore a significant shift of the 
synchrotron peak, in agreement with the rapid variation detected 
by UVOT and the drastic change of polarization angle observed 



in VLBA data simultaneously with the optical flare. The opti- 
cal/UV spectra presented here clearly show the diff'erent contri- 
bution of the thermal and non-thermal components in the various 
phases of the y-ray activity of the source, with consequent theo- 
retical implications on the modeling of the broad band spectrum 
to investigate. It will be important to continue the monitoring of 
the source from IR to UV to know if the shift of the synchrotron 
peak observed in March 2009 is a rare event or not. 

The Swift/XRT observations show no clear correlation be- 
tween the X-ray and y-ray emission. This could be due to the 
fact that X-rays and y-rays are originated by EC emission of dif- 
ferent parts of the electron distribution with different variability. 
On the other hand, a hint of harder- when-brighter behaviour was 
detected from the XRT observations suggesting the presence of 
a second component in the soft X-ray part of the spectrum that 
could be associated with the soft X-ray excess rather than to the 
variation of the SSC contribution in different activity states of 
the source. 

Moreover, the radio data show that an increase of the activity 
related to the y-ray activity started in March and peaked in April- 
May 2009, confirming that the mechanism producing the y-ray 
flare interested also the mm/cm emitting zone, with some delay 
likely due to opacity effects. 

In conclusion, the y-ray emission of PKS 1510-089 shows a 
complex correlation with the other wavelengths, leading to dif- 
ferent scenarios when modeling its broad band spectrum at dif- 
ferent epochs. The multifrequency observations presented here 
give new clues, but also offer new questions on the astrophysi- 
cal mechanisms at work in this object. Further IR observations 
can be important to test the "far dissipation" scenario in which 
the IR photons are the dominant seed photons for the IC mech- 
anism that produces the y-rays, and in this context the launch 
of Herschel (Pillbratt et al. 2010) opened a new window on the 
infrared Universe and it will be very important also for the study 
of blazars (Gonzalez-Nuevo et al. 2010). Finally, the study of 
the correlation between the y-ray activity and the variability of 
optical polarization provided new information on the jet struc- 
ture, the location and the causes of the high-energy emission in 
PKS 15 10-089. A further piece of the puzzle can be added with 
the launch of a X-ray polarimeter onboard any future X-ray mis- 
sions (e.g. the New Hard X-ray Mission, Pareschi et al. 2010) 
giving us a final answer to the real nature of the X-ray spectrum 
of this puzzling object. 
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